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ABSTRACT 

We derive the cosmic star formation history (CSFH) out to z = 1.3 using a sample of ~ 350 radio- 
selected star-forming galaxies, a far larger sample than in previous, similar studies. We attempt to 
differentiate between radio emission from AGN and star- forming galaxies, and determine an evolv- 
ing 1.4 GHz luminosity function based on these VLA-COSMOS star forming galaxies. We precisely 
measure the high- luminosity end of the star forming galaxy luminosity function (SFR > lOOMoyr -1 ; 
equivalent to ULIRGs) out to z = 1.3, finding a somewhat slower evolution than previously derived 
from mid- infrared data. We find that more stars are forming in luminous starbursts at high redshift. 
We use extrapolations based on the local radio galaxy luminosity function; assuming pure luminosity 
evolution, we derive L* oc (1 + z ) 2 - 1±0 - 2 or L* oc (1 4- z) 2 5±01 , depending on the choice of the local 
radio galaxy luminosity function. Thus, our radio-derived results independently confirm the ~ 1 order 
of magnitude decline in the CSFH since z ~ 1. 

Subject headings: galaxies: fundamental parameters - galaxies: starburst, evolution - cosmology: 
observations - radio continuum: galaxies 



1. INTRODUCTION 

Studies based on different galaxy star formation indica- 
tors (UV, optical, FIR, radio) agree that the cosmic star 
formation history (i.e. the total star formation rate per 
unit co-moving volume; CSFH hereafter) has declined by 
about an ord er of magnitud e since z ~ 1 (for a compi- 
lation see e.g. lHopkinsl T2004) . One of the major difficul- 
ties of UV/optical based tracers is the significant model- 
dependent dust-obscuration correction that needs to be 
imposed on the data. This 'dust-obscuration problem' 
may be overcome using longer wavelengths, such as the 
IR and radio regimes. However, in these cases a multi- 
wavelength approach is essential as redshift information 
and a reliable identifier of star forming (SF) galaxi es is 
required (e.g. lCaputi et al.ll2007l : [Smolcic et al.ll2008l ; S08 
hereafter). In this context the radio star formation tracer 
provides an important complementary view of the CSFH. 
First, radio emission is a dust-insensitive tracer of recent 
sta r formation (n ot affected by old stellar populations; 
see ICondonlll992l for a review). Second, interferometric 
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radio observations with ~ 1" resolution allow more reli- 
able identifications (compared to FIR and sub-mm data) 
with objects detected at other wavelengths. 

The dust-unbiased total CSFH has been constrained 
recently using MIR (24/8/mi) selected samples ob- 
tained b y deep small area surveys (CDFS, GEMS , 
GOODS: iLe Floc'h et al] 120051 iZheng et alJ l20M [20071: 
ICaputi et al.ll2007l : iBell et alJl2007T >. Small area surveys, 
however, are subject to cosmic variance. Moreover, they 
do not observe a large enough comoving volume in order 
to fairly sample rare high-luminosi ty galaxies. In this pa - 
per we use the 2D° COSMOS field (I Scoville et al]l2007al). 
and its 1.4 GHz radio observations (|Schinnerer et al.l 
2007), to derive the cosmic star formation history. In 
such a large field cosmic variance is significantly reduced 
as 2D° (1.4° x 1.4°) sample comoving volumes in the 
early universe (z ~ 1) comparable to the largest sur- 
vey in the local uni verse (SDSS - DR1; see Fig. 1 in 
iScoville et al.l l2007al ) . The physical angular size sam- 
pled by 1.4° at redshifts 0.2 - 1.1 roughly corresponds 
to a factor of 3 to 8 of the typical cluster scale length 
(~ 5 Mpc). Thus at all redshifts, such a fi eld fairly sam- 
ples re le vant structures in the un iverse (see lScoville et all 
200733 iMcCracken et al.ll2007l for a more detailed dis- 
cussions on cosmic variance in the COSMOS field). 

In the last decade several radio surveys have been 
utilized to independently derive the cosmic star for- 
mation history. However, their results are based on 
small observed areas, non-uniform rms in the final map, 
as well as a fairly non-uniform selection of SF galax- 
ies. The first derivation of th e CSFH based on ra - 
dio data has been performed by lHaarsma et alj (2000) . 
They combined three radio freque ncy observations of th e 
Hubble Deep Field (66 arcmin 2 ; iRichards et a l. 1998), 
SSA13 (7 a rcmin 2 : IWindhorst et all 1 19951). and V15 
oo arcmin , iFomalont et al.lll991t lHammer et al.lll995l) 
fields reaching 5cr radio depths of 9, 8.8, and 16 fiJy at 
the field centers, respectively. Of the total number of 
their radio-selected sources (77) only 37 were securely 
classified as star forming galaxies (based on morphology 
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and/or optical spectroscopy). Their sources reach out to 
z ~ 3, 23 have spectroscopic redshifts, and the redshifts 
for the remaining 14 h ave been estimated u sing only I 
or H and K' ba nds (see I Haarsm a et al.ll2000l for details). 
More recently, ISevmour et alJ (|2008f ) used radio obser- 
vations of the 13 H XMM-Newton/ Chandra Deep field 
(0.196D ; 4cr ~ 30 ^Jy ) to derive the CSFH. Out of 
a total of 449 radio sources they find 269 galaxies which 
they classify as star forming based on a number of crite- 
ria applicable to sub-samples of their objects (see their 
Tab. 1). About half of these galaxies have a spectroscopic 
redshift (z < 3). 

Here we utilize the 1.4 GHz VLA observations of 
the COSMOS 2D° field (VLA-COSMOS Large Project; 
ISchinnerer et al.ll2007h to overcome the above mentioned 
biases. The final mosaic has a resolution of 1.5" x 1.4" 
and a typical rms of ~ 10.5 (15) /zJy/beam in the central 
1 (2) D° making this survey to date the largest radio deep 
field at this sensitivity and angular resolution. Gi v en the 
COSMOS panchromatic data set, ISmolcic et "ail (|2008l ) 
have developed a novel method to select star forming 
and AGN galaxies based only on NUV-NIR photome- 
try. This yielded a robust statistical selection of 340 SF 
galaxies out to z = 1.3 (see Sec. 12.11) . One particular 
advantage of the VLA-COSMOS survey are the large co- 
moving volumes observed at all redshifts, thus allowing 
one to study a statistically significant sample of rare ob- 
jects, such as the most intensely star forming galaxies. In 
this paper we focus on the derivation of the CSFH, with 
emphasis on the evolution of galaxies with high star for- 
mation rates (> 100 M Q yr _1 ), using the 1.4 GHz VLA- 
COSMOS Large Project. 

We report magnitudes in the AB system, adopt a stan- 
dard concordance cosmology (H = 70, Om = 0.3, Ha = 
0.7), and define the radio synchrotron spectrum as F v oc 
v~ a , assuming a — 0.7. 

2. THE 1.4 GHZ LUMINOSITY FUNCTION FOR STAR 
FORMING GALAXIES IN VLA-COSMOS 

2.1. The star forming galaxy sample 

The sample of SF galaxies used here is presented in 
S08, and briefly summarized below. Using radio and op- 
tical data for the COSMOS field S08 have constructed 
a sample of 340 star forming galaxies with z < 1.3 out 
of the entire VLA-COSMOS catalog. The selection re- 
quired optical counterparts with i < 26, accurate pho- 
tometry (i.e. outside photometrically flagged areas), and 
a S/N > 5 detection at 20 cm, and is based on rest-frame 
optical colors (see also ISmolcic et al]|2006l ). The classi- 
fication method was well-calibrated using a large local 
sample of galaxies (~ 7, 000 SDSS "main" spectroscopic 
galaxy sample, NVSS and IRAS surveys) representative 
of the VLA-COSMOS population. It was shown that the 
method agrees remarkably well with other independent 
classification schemes based on mid-infrared colors and 
optical spectroscopy, and that, within the selection lim- 
its, it is not biased against dusty starburst galaxies. It 
is important to note that the rest-frame colors for the 
VLA-COSMOS radio sources have been shown to stay 
basically unchanged with redshift (an effect that was also 
noted by lBarger etai]l2007t iHuvnh et al"1l2008ft . 

Out of the 340 selected SF galaxies 150 have spectro- 
scopic redshifts available, while the remaining sources 
have very reliable photometric redshifts (see S08 and ref- 



erences therein). Based on Monte Carlo simulations, S08 
have shown that the photometric errors in the rest-frame 
color introduce a small, ~ 5%, number incompleteness of 
SF galaxies. Here we use the S08 sample of SF galaxies, 
statistically corrected for this effect. 



2.2. Derivation of the radio luminosity function (LF) 

We derive the radio luminosity function ($) for our 
SF galaxies in f our redshift bi ns using the standard 
l/Vrnax me thod (jSchmi dt 1968). In computing V m&yi 
we take into account both the radio and optical flux 
limits (i.e. the maximum observable redshift, z max ), as 
well as the non-uniform rms noise level in the VLA- 
COSMOS mosaic. For the latter we use the differen- 
tial visibili ty area (i.e. area! cover age, Ak, vs. rms; see 
Fig. 13 in ISchinnerer et all [2001 . For a source with 
1.4 GHz luminosity Lj its maximum volume is then 

V max (Lj) — J2k=l • Knax(Zmax> Lj). 

Several additional corrections need to be taken into 
account: i) the VLA-COSMOS detection completeness, 
ii) the fraction of sources not included in the radio-optical 
sample, and iii) the SF galaxy selection bias due to the 
rest-frame color uncertainties. 

The radio detection completen ess of the VLA - 
COSMOS survey has been derived bv lBondi etaTI (|2008l ) 
via Monte Carlo - artificial source - simulations for the 
inner ID . Artificial radio sources have been simulated 
taking into account both the flux density and angular 
size distributions. The first has been described using a 
typical broken power law, while the latter has been as- 
sumed to depend on flux density (6 mc d oc F™) 12 Com- 
parisons between the median angular sizes in different 
flux ranges and between real and simulated 1.4 GHz ra- 
dio data showed that the median angular size (O mc d) of 
radio sources is changing with flux according to a power 
law with an exponent of m = 0.5. O mc d oc F® 5 yields 
flux dependent completeness corrections. The median of 
these corrections is ~ 10% (reaching a maximum value 
of 60% in only one of the l owest flux density bins; sec 
Tab. 1 in iBondi et alJ [2008) . Assuming that their cor- 
rections, scaled for the differences in radio sensitivity, 
hold for the full 2D°, we utilize them to correct our LFs 
for the detection incompleteness (fdet)- 

To correct for objects located in photometrically- 
flagged regions in the optical images, or not identified 
with an optical counterpart, and thus not included in 
the radio-optical sample, we construct a correction curve 
as a function of total 1.4 GHz flux density which yields 
an average correction of ~ 30% (see Fig. 23 in S08; ffl ag ). 
Hence, in the i th luminosity bin the comoving space den- 
sity and its corresponding error (<7i), are computed 
by weighting the contribution of each galaxy by the two 
correction factors, fdot and ffl ag , which were obtained by 
linearly interpolating the two correction curves described 
above, respectively, at the total flux of the given, j th , 

12 The angular size dependence of radio sources on their flux 
density has been observed in numerous radio surveys, and it is 
mainly due to the radio luminosity - size relation. At fainter flux 
densities, intrinsically fainter (and therefore smaller) radio sources 
are preferentially sampled. 
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The selection bias due to the rest-frame color uncer- 
tainties is accounted for via Monte Carlo simulations. 
As described in S08, the rest-frame color error distribu- 
tion is simulated using a Gaussian function with a = 0.1. 
These errors are then added to the rest-frame color ob- 
tained from the SED fitting. The SF galaxies are then 
selected, and the LF is derived as described above. This 
process is repeated 300 times, hence we obtain 300 real- 
izations of for each luminosity bin. We take the 
median values as representative. 

2.3. The radio luminosity function 

The LFs for our SF galaxies for the 4 chosen redshift 
bins are shown in Fig.[TJ and tabulated in Tab.Q] In each 
panel in Fig. Q] we show the analytical form of the locally 
derive d 20 cm LF for SF galaxies given by I Sadler et al.1 
(2002, see also Sec. 12.41) . It is worth noting that although 
the 2D° COSMOS field samples a relatively small comov- 
ing volume at the lowest redshifts and only a photometric 
identification of SF galaxies has been used, our LF in the 
lowest redshift bin agrees remarkably well with the local 
LFs that were derived using all-sky radio surveys (NVSS) 
combined with good quality optical spectroscopic data 
(SDSS, 2dF) to identify SF galaxies. For comparison, 
in the first redshift bin we also indicate other local radio 
LFs t hat exist in the literature (|Condonl [l989: Bes t et all 
|2005() . We discuss the implications of the differences be- 
tween the various local radio LFs on the star formation 
history further below (Sec. I2.4.T]) . 

In all panels we show the volume densities of 20 cm 
radio sources derived by lHaarsma et al.1 (|2000t gray as- 
terisk), corrected to the current cosmology. Haarsma 
et al. used 37 star forming galaxies to derive these LFs 
(38% of these had approximate redshifts derived from I- 
or HK'- band magnitudes, the others had spectroscopic 
redshifts) . Their data points in each redshift range agree 
fairly well within the error-bars with our results. Note 
that due to the almost one order of magnitude larger 
sample of star forming galaxies used here, the error-bars 
of the VLA-COSMOS LFs are significantly smaller. 

In the higher redshift panels w e also compare our LFs 
with the total IR LFs derived bv lLe Floc'h et all (|2005l . 
hereafter LF05) based on a 24 /im selected sample in the 
CDFS (Chandra Deep Field South; top right and both 
bottom panels in Fig. [l| . The total IR luminosity was 
converted to 1. 4 GHz lum inosity using the total IR - ra- 
dio correlation (Bell 2003) , which has an intrinsic scatter 
of ~ 0.26 dex, shown by horizontal error bars in Fig. [TJ 
The IR LFs were re-scaled to our redshift ranges either by 
combining two narrower redshift bins given in LF05 or by 
scaling a given comoving density using the evolution pa- 
rameters, and their corresponding errors, given in LF05. 
There is an excellent agreement between the 1.4 GHz and 
IR LFs. Note also that the VLA-COSMOS LFs constrain 
much better than the IR data the high-luminosity end, 
i.e. the most intensely star-forming galaxies. Interest- 
ingly, the VLA-COSMOS star forming galaxies in our 
two highest redshift bins (z > 0.6) show an extended 
high-luminosity (h 1A > 2 x 10 24 W Hz" 1 ) tail. We 
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Fig. 1. — 1.4 GHz luminosity functions (LFs) for star form- 
ing galaxies in the VLA-COSMOS survey, shown for four redshift 
ranges (filled blue squares) are presented in each panel. The num- 
ber of galaxies in each redshift bin, statistically corrected for se- 
lection uncertainties (see text for details), is also indicated in each 
panel. The local 20 cm LFs for star forming galaxies are shown 
in th e top -left panel {C ondon 1989; S adler et atll2002l ; IBest et al.1 
[20051; the ISadler et al.l 120021 LF — solid orange curve — is shown 
in all panels). In all pan els we show the volume densities derived 
by lHaa rsma et al. (200(3, g rav asterisk), corrected to the current 
cosmology. In the top right, and bottom panels, the total IR LFs 
dLe Floc'h et al. 2005) for the corresponding redshift ranges (open 
circles) is also shown (see text for details). The total IR luminosity 
was conver ted to 1.4 CHz luminosity using the correlation given in 
iBeTH lr2003h . 

cannot exclude some possible contamination from AGN, 
which are much more numerous than SF galaxies at these 
high radio luminosities (see Fig. 17 in S08). However, a 
similar excess of SF galaxy volume densities at the high 
luminosity end has been recently found by Cowie et al. 
(2004), who have used the HDF-N and SSA 13 fields 
to derive the radio LF for spectroscopically identified SF 
galaxies (106 SF galaxies, z < 1.4). 

2.4. Towards the derivation of the cosmic star 
formation history 

The comparison of our derived LFs (see Fig.[T|) with the 
local 20 cm LF shows a strong evolution with look-back 
time. The evolution of a given population of objects is 
usually parameterized with monotonic density and lumi- 
nosity evolution. However, as the VLA-COSMOS data 
do not allow the derivation of the LF out to, and fainter 
than the characteristic luminosity (i.e. the 'knee' of the 
LF), a full determination of the evolution is not possible 
with these data, and we choose to parameterize it using 
pure luminosity evolution. However, as we show below, 
its quantization depends fairly strongly on the choice of 
the local LF, as those presented in the literature are not 
exactly the same. 

2.4.1. The local 20 cm luminosity functions 
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The two commonly used analytical for ms for th e loca l 
20 c m radio LFs h ave b een presented in ICondonl (|1989l ) 
and ISadler et all (|2002l) . Cond on et al. use a hyp er- 
bolic parameterization (see also ICondon et al.1 l2002f ) of 
the form: 



log $i.4GHz = 28.83 - - log Ll.4GHz + Y 



B 



(logLi,4GHz ~ X) 
W 2 



1/2 



(2) 



where Y = 2.88, B = 1.5, X = 22.108, W = 0.667 
(corrected to the current cosmol ogy and the base of 
dlogL, compared to d mag given in lCondonlll989f h These 
parameters have been derived using a sample of 307 
spiral and irregular g alaxies from the Rey ised Shapley- 
Ames Cat alog (RSA: ISandage et al.lll98lT ) observed at 
1.49 GHz (|Condonlll987l see Fig. 12^ 

Local 20 cm luminosity functions 




RSA galaxy sample 
• 2dFGRS + NVSS galaxy sample 

- Condon 1989 (RSA) 
— Sadler et al. 2002 (RSA+2dF) 
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Fig. 2. — Illustration of the differences in the shapes of the 
local 1. 4 GHz radio luminosity function s (LFs) found in the lit- 
erature (Condon 1989; Sadler et al. 2002, corrected to the current 
cosmol ogy) . A l so sho wn are the data points used for the analytical 
fits in Condon (1989, Revised-Shapley-Ames Catalogue; RSA) and 
ISadler et al.l ( 120021 . NVSS/2dFGRS and RSA samples). Note that 
the major difference in the representation of the local radio LFs 
is the analytical form used to fit the data (see eqs. fS] & |3J that 
particularly affects the position of the LF break (i.e. 'knee'). 



The LF given in ISadler et al.l (|2002l ) takes on the form 
of a combined power-law and Gaussian distribution given 
by the following a nalytic function (first proposed by 
ISandage et ai1ll979l ): 



$(L) = $* 



L 



exp 



1 



log(l + 



L_, 



(3) 



with a = 0.84, a = 0.94, = 22.9 x 10" 3 Mpc" 3 , 
and L* = 1.95 x 10 19 W Hz" 1 (scaled to the cosmol- 
ogy u se d here, and t o the base of cilogL; see iHopkinsI 
120041) . ISadler et alj (|2002t ) have used 204 SF galaxies 



drawn from the Two Degree Field Galaxy Redshift Sur- 
vey (2dFGRS) and the 1.4 GHz NRAO VLA Sky Survey 
(NVSS) to derive their LF (see Fig. [3]). However, to ob- 
tain the best fit parameters to eq. [3l they combined this 
sample (that constrains well the high luminosity end of 
the LF; see Fig. [2]) with the RSA galaxy sample (in order 
to sample the low luminosity end of the LF). 

The differences between the two local LFs are illus- 
trated in Fig. [21 There is a discrepancy between the 
two analytical representations of the radio LFs at both 
the high and low luminosity ends. However, the dis- 
crepancy affecting most severely the star formation rate 
density, that we aim to derive here, is the different po- 
sition of the LF's turn-over point (L*) given by the two 
analytic forms. This yields a difference of 10-50% in the 
star formation rate density integral (see e.g. Fig. [3] and 
Fig. [5]) as the luminosity range encompassing the turn- 
over point (10 19 < Li. 4G Hz < 10 23 W Hz -1 ) contributes 
the most (~ 95%) to the integral. It is impo rtant to note 
that t he 2dFGRS and NVSS data used by ISadler et aD 
(2002) sample more precisely the high-luminosity end 
(see Fig. [2]) o f the local radi o LF when compared to the 
RSA sample ()Condonlll989h . Hence, in Sec. GO we will 
use the former one to constrain the evolution of our most 
intensely SF galaxies (i.e. L 1A > 2.34 x 10 23 W Hz" 1 ). 
As the total star formation rate density in each redshift 
range is derived by integrating under the evolved LF 
curve (multiplied with star formation rate; see below) 
down to the faintest 1.4 GHz luminosities, we will take 
this difference of the local 20 cm LFs into account in 
further analysis. 

2.5. The evolution of star forming galaxies 

We parameterize the evolution of the VLA-COSMOS 
SF galaxy LF by pure luminosity evolution: 



$ Z (L) = $, =0 



L 



(4) 



where oll is the characteristic luminosity evolution pa- 
rameter, and $> z is the luminosity function at redshift 
z. We derive the evolution by summing the x 2 distri- 
butions obtained for a large range of fixed «l in each 
redshift bin (excluding our first - local - redshift bin). 
The uncertainty in aj, is then taken to be the 2a error 
obtained from the x 2 statistics. Our results yield a pure 
luminosity evoluti on with i) aj, = 2.1 ± 0.2, when the 
ISadler et all (I2002D local LF is used, and ii) a L = 2.5±0.1 
when the ICondonl (|1989l ) local LF is used. The different 
evolution parameters are a natural consequence of the 
different slopes of the two local LFs in the luminosity 
range that is constrained by the VLA-COSMOS data 
(s ee Figs. El and II 

lHaarsma et al.l ()2000D have found that a pure lumi- 
nosity evolution with qj, w 2.74 is a good representa- 
tion of the evolution of their radio-selected SF galax- 
ies [no uncerta inties were asso ciated with this estimate] . 
They used the ICondonl (|1989l ) local LF as the basis for 
deriving their evolution. Given i) their much smaller 
sample size, and ii) that their LF, when corrected for 
cosmology, agrees well with the one derived here (see 
Fig. [J ), these two re sults are in good agreement. Further, 
ICowie et al.l (|2004|) find an e volution of the SF galaxy 
LF (using the ICondonl 1 19891 local LF) consistent with 
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qj, of 3. The SF LFs, as well as the evolu tion, derived 
here agree with their findings (c.f. Fig 3 in ICowie et al.1 
2004). Our results are also consistent with th e overall 
evolut ion of star forming galaxies obtained by [Hopkins 
(|2004 a L = 2.7 ± 0.6), that has been shown to fit well 
the evolution of radio selected star form ing galaxies at 
low redshifts (z < 0.5: lAfonso et al.ll2005l Phoenix Deep 
Survey). 

3. THE COSMIC STAR FORMATION HISTORY (CSFH) 
3.1. The total cosmic star formation history 
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Fig. 3. — Luminosity density for VLA-COSMOS star forming 
galaxies (filled squares) in 4 redshift bins. The solid and dashed 
curves correspon d to the best fit LFs in ea ch redshi f t bin using the 
local LF given in S adler et~aTl [|2002h and ICondoiil (fT989h ■ respec- 
tively. Note that the first appears to describe the bright-luminosity 
end more consistently with our data. 

As the star formation rate density (SFRD) is derived 
by integrating the luminosity density, in Fig. [3] we show 
the luminosity density for our 4 redshift bin s. The 
shown curves are the two local 20 cm LFs (jCondonl 
1989; Sa dler et al.ll2002f) . purely evolved in luminosity, 
and best fit to the VLA-COSMOS data in each redshift 
range. Prior to integration, we convert the 1.4 GHz ra- 
dio luminosity to star formati on rates (tb , in M^yr -1 ), 
using the calibration given in [Belli (|2003l) . based on the 
total IR - radio correlation: 



tf> [Mq yr" 



5.52 x 10~ 22 L 

5.52x10 



0.1+0.9(Li. 4 /L c ) 



1.4 , 



Li.4 > L c 
L1.4 < L c 



(5) 

where L c = 6.4 x 10 21 W Hz^ 1 , and L 1A is the 1.4 GHz 
radio luminosity in units of W Hz -1 . This calibration 
uses a Salpeter initial mass function (IMF oc M -2 35 ) 



from 0.1 - 100 M a . After the conversion, we compute the 
star formation rate density (SFRD) for a given redshift 
bin as J ip (L) <E> Z (L) dL, where $ z is the evolved radio 
LF best fit to the data in each redshift range (see curves 
in Fig. El). 

In Fig. |4] we first show the SFRD, obtained by numeri- 
cally integrating the VLA-COSMOS data (open squares) 
within the sampled luminosity range, with no attempt 
made to extrapolate towards the faint or bright lumi- 
nosity ends using the evolved local LF. This eliminates 
any assumption, and yields robust lower limits purely 
obtained from the data. The dependence of the derived 
SFRD on the faint /bright end extrapolation is illustrated 
by the (solid and dashed) curves in Fig. [H which were ob- 
tained using the average pure luminosity evolution of the 
two local LFs derived in Sec. I2.51 and integrating over the 
entire SFRD curve. Hence, the largest uncertainty in the 
derivation of the SFRD based on radio data arises from 
the uncertainty in the shape of the local radio LF (as- 
sumed not to change with redshift), and the associated 
extrapolation below the faintest luminosity sampled by 
the data. 



0.100 




d, 0.010 



0.001 



-i— i— i— |— i— i— i— r-i— i— i— r 

Local 20cm results: 

• Sadler et al. (2002) 

• Condon et al. (2002) 
X Machalski & Godlowsk! (2000) 
© Serjeant et al. (2002) 
V Mauch & Sadler (2007) 
O Condon (1989) ^ ' ' 



ft* 



VLA-COSMOS lower limits □ 




0.0 



0.2 



0.4 



0.6 



1.0 



1.2 



Fig. 4. — Star formation rate density (SFRD) as a function 
of redshift. Open squares represent the SFRD derived from pure 
VLA-COSMOS data, without any extrapolations towards the faint 
or bright luminosity end. Therefore, they present strict lower lim- 
its, and are indicated with arrows. The two curves shown corre- 
spond to the average evolution of the VLA-COS MOS SF galax- 
ies derived using two different local radio LFs llCondonl 119891 ; 
ISadler et all 12002) . Local SFRDs derived by various authors are 
also shown, and indicated in the panel. 



We compare our SFRD results (given in Tab. [2]), ob- 
tained by integrating over the best fit evolved local LF 
in each redshift range (see curves in Fig. [3|) , with other 
radio-based estimates in Fig. [5] The expected steep de- 
cline in the star formation rate density since z ~ 1 is re- 
produced by our data. Our result s are consistent within 
the errors wi th the results from lHaarsma et al.l (|200Ct 
who used the ICond"onlll989l local LF), although their re- 
sults are on average higher. Note also that our statistical 
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TABLE 1 

1.4 GHZ RADIO LUMINOSITY FUNCTIONS FOR 
VLA-COSMOS STAR FORMING GALAXIES 

(H = 70, U M = 0.3, n A = 0.7). 



r— 

redshift 1 


Ll.4GHz 






range 


[WHi 


J 


[Mpc~ 3 dex _1 ] 


| 


1.0 ■ 


1Q 22 


7.0 ± 1.5 • 


10" 4 


U.l < 2 n U.ou 


2.0 ■ 


10 22 


2.9 ± 0.7 ■ 


10" 4 




4.0 ■ 


10 22 


3.5 ±0.6 • 


io~ 4 




8.9 ■ 


10 22 


5.0 ±2.1 • 


10 -5 




1.8 ■ 


10 23 


3.2 ±1.6 • 


10 -5 




3.2 ■ 


10 23 


1.4 ± 1.0 • 


10- s 




3.5 ■ 




2.4 ±0.9 ■ 


10~ 4 


U.oD < ^ S U.O 


6.3 ■ 


10 22 


1.5 ± 0.4 • 


10" 4 




1.0 ■ 


10 23 


1.3 ± 0.3 • 


10~ 4 




1.6 ■ 


10 23 


5.8 ±1.6 • 


10~ 5 




2.2 ■ 


10 23 


9.0 ±6.4 • 


10~ 6 




3.5 ■ 


10 23 


1.5 ±0.7 • 


10~ 6 




1.4 ■ 


10 M 


1.2 ±0.2 • 


io- a 


0.6 < z < 0.9 


2.8 ■ 


10 23 


4.8 ±0.9 • 


10 -5 




5.6 ■ 


10 23 


1.7 ±0.5 ■ 


10 -5 




1.1 ■ 


10 24 


5.8 ±2.6 • 


10~ 6 




2.5 ■ 


10 24 


5.2 ±2.3 • 


10" 6 




5.0 ■ 


10 24 


2.0 ±1.4 • 


10~ 6 




3.2 ■ 


W 2i 


9.2 ±1.5 • 


10 -t> 


0.9 < z < 1.3 


6.3 ■ 


10 23 


2.4 ±0.5 • 


10 -5 




1.3 ■ 


10 24 


1.8 ±0.4 • 


10" 5 




2.8 ■ 


10 24 


3.2 ±1.4 • 


10~ 6 




5.6 ■ 


10 24 


1.2 ±0.8 • 


10" 6 




1.0 ■ 


10 25 


2.0 ±1.0 • 


10" 6 




1.8 ■ 


10 25 


1.0 ±0.7- 


10~ 6 



TABLE 2 

Star formation rate density (d erived using the 
calibration given in ibell][2003t see eq. [5} from 
1.4 GHz VLA-COSMOS data 
(H = 70, Q A/ = 0.3, U A = 0.7). 



redshift I SFRD [ M^yr" 1 ] 

range ICondonl jl989f) LF ISadler et all ll2002f) LF 



CSFH uncertainties (bold red cro sses in Fig. are si. 
nificantly smaller compared to the lHaarsma et aT] (200 1 
results, as a result of the larger sample utilized here. 
Our results are also qual i tative ly consistent with those 
obtained by llvison et al.l (|2007f ) , based on radio image 
st acking of UV se le cted g alaxies in the AEGIS20 survey. 

ISevmour et~aTl (|2008h have used VL A/MERLIN 
radio frequency observations of the 13 H XMM- 
Newton/ Chandra Deep field to derive the cosmic star 
formation history. Thei r findings are qualitat ively con- 
sistent with those from lHaarsma et al.1 (l2000l. sho wn in 
Fig. [5]) when they use the iMauch &: Sadlerl (|2007f) local 
LF evolved in luminosity with an a-priori set value of 
Qi = 2.5 (note that this loc al LF has a lower normaliza- 
tion than lSadler et al.ll2002t see the local results in Fig. [4] 
and Fig. [5]). Further, they used a different radio lumi- 
nos ity to SFR relation that is consistent with 0.84 times 
the | Belll (|2003h calibration used here (see ISevmour et al.l 
120081 for details). Hence, this implies that their results 
are significantly higher than the VLA-COSMOS results. 



Further, given the differenc es between the local radio 
LFs outlined in S ec. [2XT1 i f ISevmour et all (|2008| ) had 
chosen to use the Condon (1989) local LF with other- 
wise the same assumptions, they would have obtained 
even higher SFRD values. We believe that the main rea- 
son for the diff e rence s between our results and those by 
ISevmour et al.l (|2008l ) is likely a combination of i) their 
a-priori assumed evolution of the local LF, contrary to 
constraining the evolution by their data and ii) the inclu- 
sion of a fraction of lower power (i.e. radio quiet) AGN in 
their star forming galaxy sample, while our sample may 
over-subtract composite objects (see also Sec. 0}. 

In Fig. M we compare the VLA-COSMOS derived 
CSFH data with results from previous studies based on 
a range of SF estimators - UV, optical, FIR, total IR, 
and radio. A luminosity-dependent obscuration co rrec- 
tion was used where necessary fsee iHopkinsI 12004 and 
references therein). Overall, our derived CSFH agrees 
with the general trend of a rapid decline by almost an 
order of magnitude in the cosmic star formation rate den- 
sity since z ~ 1 . A possibly s lower decline is suggested by 
our data if the lSadler et al.l (|2002f ) LF is used. However, 
as the uncertainties due to the faint-end extrapolation 
are significant, no robust conclusions can be made at 
this point. 
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Fig. 5. — Comparison of radio derived cosmic star formation 
histories. The different samples are indicated in the pa nel. The 
bold r ed circles show the VLA-COSMOS results, when the lCondonl 
(1989) local LF is u sed; open red circles show the results when the 
ISadler et al. ( 2002) local LF is used. Results based on other radio- 
based studies are also shown, and indicated in the panel. 



3.2. The CSFH of intensely star forming galaxies 

The VLA-COSMOS SF galaxy sample constrains well 
the high end of the LF for SF galaxies. Given the 2D° 
VLA-COSMOS field the comoving volume sampled up to 
z = 1.3 is - 1.3xl0 7 Mpc 3 , corresponding roughly to the 
volume observed locally by SDSS (DR1). Thus, for the 
first time this allows a robust derivation of the CSFH 
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for galaxies forming stars at rates of > 100 M© yr _1 
out to z = 1.3. Such radio selected galaxies are equiva- 
lents to the ultra-luminous IR galaxies (ULIRGs, Ltr > 
10 12 L ), and it is noteworthy that the VLA-COSMOS 
survey is sensitive to a complete sample of these galaxies 
out to z ~ 1 (see Fig. 16 in S08). 

In order to derive the evolution of the SFRD at the 
high-luminosity end, we integrate the SFRD curve, ob- 
tained from the best fit pure radio luminosity evolution 
in each redshift bin (see Fig. [3]) , only for our SF galaxies 
that have L1.4 > 2.34 x 10 23 W Hz -1 , which corresponds 
to Lir > 10 12 Lp) g iven the adopted total IR - radio cor- 
relation (p 3c 1 2003). For this we use the local LF given by 
I Sadler et alj p002) as it appears to be b etter sui t ed for 
the high-luminosity end compared to the ICondonl (|1989f ) 
LF (see Fig. [5]). Note that for these highly luminous 
galaxies the extrapolation uncertainties are not as signif- 
icant as for the overall SFRD, as this sample is almost 
complete in all three high redshift ranges. A small ex- 
trapolation to the faint end, given the form of the evolved 
local LF, is necessary only in the last redshift bin. For a 
consistent comparison between our radio and IR (LF05) 
results we convert the total IR luminosity to star fo rma- 
tion ra tes consistently using the calibration given bv lBelll 
(2003). The results are shown in Fig. [7] The evolution 
of our star forming ULIRGs is consistent with the lower 
envelope derived by LF05. However, it is marginally flat- 
ter, suggesting a slower evolution of star forming ULIRG 
galaxies since z — 1.3. This will be further discussed in 
the next section. 



4. DISCUSSION 

4.1. The evolution of the most intensely star forming 

galaxies 

Making use of our large statistical sample of radio- 
selected star forming ULIRGs complete out to z ~ 1 
we have derived the CSFH of the most intensely star 
forming galaxies (> 100 M Q yr ) out to z = 1.3. 
Our evolution of the cosmic star formation rate in star 
forming ULIRGs qualitatively agrees with previous MIR- 
based results (LF05; see Fig. [7]). However, we find a 
slightly slower evolution than predicted by the MIR re- 
sults. The major cause for this difference is currently 
unclear, nonetheless there are likely three effects that 
may contribute: 

i) No attempt has been made to minimize the AGN 
contamination in the 24 fmx- selected sample possibly 
causing an overestimate in the MI R-derived SFRD evo - 
lution for ULIRGs. For example, iCaputi etlrtl (|2007f ) 
have found > 10% of 24 ^m-AGN at z ~ 1, and a 
factor of 2 more at z ~ 2, suggesting that the AGN 
fraction in MIR selected samples increases with redshift. 
The AGN fraction in MIR samples may also be a func- 
tion of stellar mass. Although at higher redshift than 
analyzed he r e (im plying a different cosmological era), 
iDaddi et~a l. (2007) have demonstrated that at z ~ 2 the 
MIR AGN fraction is indeed a function of stellar mass, 
and reaches ~ 50 - 60% for masses > 4 x 10 10 M . The 
median stellar mass of the VLA-COSMOS star forming 
galaxies (obtained via SED fitting; see S08 for details) is 
~7x 10 io M . 

ii) Particular care was taken to separate the VLA- 
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Fig. 6.— Comparison of the VLA-COSMOS derived cos- 
mic star formation history (red symbols; see Fig. \E\ for details) 
with other-wavelength based results. The compilation of CSFHs 
contains UV-, Ha—, FIR-, and X-ray- based results, corrected 
for dust-obscuratio n using lumino sity-dependent corrections (open 
gray squares; see Hopkins 2004] f or details). Thick gray dia- 
monds denote local radio estimat es (Machalski & Godlowskjl l2000l ; 
ICondon et "all 12003 : ISadler et aTll20021; ISerieant et al.1120021) . For 
all radio data th e 1.4 GHz l uminosity to star formation rate cal- 
ibration given in Bell (2003, sec cq. is used. The gray-shaded 
area shows the cosmic star formation history derived by LF05, 
based on 24 (im data. The open gray circles denote i) the Ho de- 
rived SFRD at z = 0.24 in the COSMOS field, corrected for both 
dust-obscuration and AGN contribution (Shiova et al. 200<|), and 
ii) the dust-obscuration corrected SFRD at 2 = 1.2 derived us- 
ing [OH] emission gal axies in the Subaru Deep Field (SDF; see for 
Takahashi et al. 2007 details). 



COSMOS population into SF and AGN galaxies. 
Nonetheless, it has to be noted that some uncertainty, 
due to AGN contamination as well as incompleteness of 
the star forming galaxy sample exists, especially on the 
galaxy-b y-galaxy and compo site (SF plus AGN) galaxy 
level (see lSmolcic et al.ll2008l for details). It is also worth 
noting that two significant large scale s tructure com- 
ponen ts exist in our highest redshift bin ([Scoville et al.l 
l2007al Scoville et al., in prep) that may affect the frac- 
tion of SF and AGN galaxies present in this particular 
redshift range. 

iii) The local IRAS total IR and the Sadler et al. radio 
LFs, used for the derivation of these results, are not per- 
fectly similar. There may also be a volume density excess 
of SF galaxies with high radio luminosities at z > 0.6, 
compared to the Sadler et al. LF evolved only in lumi- 
nosity. In addition, the total IR - radio correlation, 
used to select ULIRGs from our radio sample, carries 
its own uncertainty and intrin sic astrophysical scatter 
(jBressan et al.ll2002t lBell2003l ). Therefore, it is not im- 
mediately obvious whether the same results would be ex- 
pected based on both - radio and MIR - star formation 
indicators. 

In summary, it is encouraging that the same qualita- 
tive behavior of the evolution of the ULIRG population 
is observed with the two independent, radio and MIR, 
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Cosmic stor formation history (ULIRGs) 




Fig. 7. — Cosmic star formation history (CSFH) derived from 
VLA-COSMOS data for the overall population (red symbols as 
in Fig- EJ, and o nly for star forming ; ULIRGs (magenta symbols; 
derived using the Sadler et al. 2002 local LF evolved in luminos- 
ity in each redshift bin). Also shown is the CSFH for the entire 
(light-gray shaded curve), and the ULIRG population (dark-gray 
shaded curve), derived from the evolved total IR luminosity func- 
tion (LF05). The IR- and radio- based star formation rates have 
been put on the same relative scale (see text for details) . Note that 
the VLA-COSMOS data suggest a slightly slower evolution of the 
star forming ULIRG population with redshift. 

SFR indicators. However, further dedicated studies of 
the details will prove most interesting in understanding 
the quantitative differences seen in Fig. 

4.2. The rapid decline of the CSFH since z ~ 1 

Our overall CSFH (Fig. agrees well with past find- 
ings, when these are corrected for dust-obscuration as 
needed. This verifies the assumptions about large dust- 
obscuration corrections required, especially for short- 
wavelength (e.g. UV) star formation tracers. Our radio 
data independently confirm the ~ 1 order of magnitude 
decline in the cosmic star formation rate since z ~ 1. 

Based on UV and IR based SF R/morphology studies 
(IWolf et al.l(200llBell et al.|[200l iMelbourne et alJ l2005t 
lHammer et al.ll2005t iZamoiski et al.ll2007T> this rapid de- 
cline in the overall cosmic star formation history is ex- 
pected to_be_ohiyeii_by normal spiral galaxies. For ex- 
ample, iBell et a l. (2005) have performed a detailed mor- 
phological study of a galaxy sample (z = 0.7) with SFRs 
> 10 M Q yr _1 . They have demonstrated that physical 
processes that do not substantially affect galaxy mor- 
phology, such as minor mergers, gas consumption and 
weak interactions with satellite galaxies, may be most 
importan t for the rapid decline in the overall CSFH (see 
also e.g. iSomerville et al.ll2001l ). 

The sample analyzed here, however, is most sensitive 
to ULIRGs (SFR > 100 Moyr" 1 ), and those are the 
systems a priori expected to be a reflection of galaxy 
merging (based on loc al ULIRG morphology studies; 
ISanders fc Mirabel Il996| ). This implies that the rapid 
decline in CSFH of our ULIRGs (Fig. [7]) may be more af- 
fected by galaxy mergers than the overall CSFH decline 



since z ~ 1. Intriguingly, our derived LF evolution is 
very similar to the evolution of the gala xy close-pair frac- 
tion derived for bright galaxies (> L v ; iKartaltepe et al.1 
I2007T) in the COSMOS field, well described with a power 
law with an index of 3.1 ± 0.1 (or 2.2 ±0.1 wh en pure lu- 
mino sity evolution is considered; see also e.g. lLotz et all 
2008). This strong evolution of the bright close-pair frac- 
tion, combined with a similarly strong evolution of the 
CSFH of ULIRGs derived here suggests that major merg- 
ers may play an important role in the rapid decline of 
the CSFH since z ~ 1 for the most intensely star form- 
ing galaxies. Thus, the observed steep evolution of our 
ULIRG population may be a good record of merger rate 
evolution, combined with gas content evolution. 



5. SUMMARY 

We have derived the cosmic star formation history 
out to z = 1.3 using to date the largest sample of 
radio-selected star forming galaxies observed at 1.4 GHz 
(20 cm) in the VLA-COSMOS survey. The large in- 
crease in the number of radio selected SF galaxies out 
to high redshift, compared to previous studies, allowed 
us to constrain well the evolution of the 1.4 GHz lumi- 
nosity function for radio-selected star forming galaxies, 
as well as to reduce significantly the statistical uncer- 
tainties of the radio-derived CSFH. We find that the un- 
certainties are ruled by the differences in the shape of 
the local radio LFs present in the literature. A pure 
radio luminosity evolution of VLA-COSMOS star form- 
ing galaxies is wel l described with L* a (1 + z) 21±0 2 , 
when evolving the I Sadler et al.l (|2002l) local LF, or w ith 
cx (1 + z ) 2 - 5±01 when evolving the lCondonl (|1989f ) lo- 
cal LF. Although encompassing a relatively broad range, 
both values are consistent with pre viously derived evo- 
lution of star forming galaxies (e.g. ICondon et all 120021 : 
lHookinsll20r)i iCowie et al.ll20t)l lAfonso et al.ll2005l )~ 

Our overall CSFH agrees well with past findings, 
when these are corrected for dust-obscuration where 
needed. This verifies the assumptions about large dust- 
obscuration corrections required, especially for short- 
wavelength (e.g. UV) star formation tracers. Making 
use of our large statistical sample of radio-selected star 
forming ULIRGs complete out to z ~ 1 we have robustly 
constrained the high-end of the SF galaxy LF at different 
cosmic times. Using these we have derived the CSFH of 
the most intensely star forming galaxies (> 100 M Q yr , 
i.e. star forming ULIRGs) out to z — 1.3. We find an, 
on average, slower evolution of the cosmic star formation 
rate in star forming ULIRGs than predicted by MIR re- 
sults consistent with the fraction of star forming galaxies 
in MIR samples likely becoming lower with increasing 
redshift and/or stellar mass. 
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